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1 Completed Work 

Essentially all the work proposed in the proposal has been completed or adequately 
addressed. Additional work on heat-exchanger aerodynamics not in the original pro- 
posal was also conceived and executed with the agreement of the Technical Monitor. 

A summary of the specific tasks accomplished is given below. The corresponding 
detailed results have been transmitted to Dryden personnel via e-mail as soon as 
they became available. All e-mail transmissions have been saved, and the key ones 
are given in the Appendices, along with relevant plots, figures, and papers which were 
generated. 


1.1 Airfoil design 

The Apex- 16 airfoil was designed specifically for the APEX test vehicle. It is in- 
tended to be representative of airfoils required for lighweight aircraft operating at 
extreme altitudes, which is the primary research objective for the APEX program. 
In the course of designing this airfoil, Fairly extensive studies were made over the 
Mach number and Reynolds number ranges of interest. Also considered were the 
airfoil thickness, pitching moment, and off-design behavior. Limited use was made 
of the optimization driver LINDOP [1] to fine-tune the final design and resolve the 
numerous conflicting requirements. Appendix B shows the Apex-16 airfoil geometry, 
coordinates, and computed performance polars. 


1.2 Study of airfoil constraints on pullout maneuver 

The maximum ceiling parameter M 2 Cl value achievable by the Apex-16 airfoil (or any 
airfoil for that matter) was found to be a strong constraint on the pullout maneuver. 
It was concluded that if data is to be acquired in level flight, then some sort of lift 
augmentation would be required, since any airfoil has little or no excess-lift capability 
at its ceiling condition (by definition). A workable alternative approach which was 
identified is to use wing lift to achieve pullout, but then acquire data in windup turns 
to achieve the target Cl- 


1.3 Selection of tail airfoils 

Several candidate airfoils for the tail were examined. The primary goal was pre- 
dictable behavior at low Reynolds numbers, and good tolerance to flap deflections. 
The NACA 1410 and 2410 airfoils (inverted) were identified as good candidates. 


1.4 Examination of wing twist 

A number of issues related to wing twist were examined. These included: simplicity, 
obtaining a uniform local c t across the test section, obtaining a high overall C'i max 
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for the pullout, avoiding tip stall. Both washin and washout options were examined 
computationally. In the end it was decided that a simple flat wing was a reasonable 
compromise between all the requirements, and allowed the use of one mold for both 
wing halves — a significant cost saving. 


1.5 Test section instrumentation and layout 

It was decided that the instrumentation for the test section would consist of surface 
pressure taps, wake rakes, surface-mounted microphones, and skin-friction gauges. 
Using multiple wake rakes was deemed desirable to verify spanwise uniformity, but 
was found to be too demanding of the limited data bandwidth available. Using a 
single wake rake mounted on a mechanical traverse was rejected due to weight and 
complexity. It was decided to use a single fixed wake rake for accurate measurement of 
the wake momentum defect, and multiple integrating rakes to verify spanwise unifor- 
mity. Several test section layouts were designed, with varying degrees of compactness. 
One of these was selected by Dry den personnel for use on APEX. 

The wake rake was designed specifically to capture the wakes anticipated at test 
altitudes. A design guide for the integrating rakes was also prepared. The layouts 
and design guide are attached as Appendices B and C. 

1.6 Integrated airfoil/heat-exchanger tests 

A modest wind tunnel test was performed for an integrated airfoil/heat-exchanger 
configuration, which is currently on Aurora’s Theseus aircraft. Although not directly 
related to the APEX tests, the aerodynamics of heat exchangers has been identified as 
a crucial aspect of designing high-altitude aircraft, and hence is relevant to the ERAST 
program. The computational studies and tunnel tests have proven the validity of 
integrating the heat exchanger with the wing airfoil to achieve surprisingly low drag 
levels at low Reynolds numbers. The results have already appeared in the Jan-Feb 
’96 issue of the Journal of Aircraft [2]. This paper is attached as Appendix D. 
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Figure 1: Integrating rake layout and dimensions 


1 Basic Relations 


An integrating rake consists of an array of pitot tubes all feeding into a common reservoir, 
as shown in Figure 1. Assuming that the static pressure across the entire wake is constant 
and equal to the edge value p e , the wake velocity profile u{y) will produce a y-dependent 
total pressure seen by the tubes: 

po(y) = \pu 2 {y) + Pe (!) 

= \pu 2 (y) + p 0oo - \pu\ ( 2 ) 


The non-uniform p 0 {y) will produce a mass flow riii in the i’th tube. Assuming fully- 
developed Poiseuille flow in each tube, this mass flow is proportional to the driving pressure 
difference p„ - p r , where p r is the reservoir pressure to be measured. 


m,- 


7T d 4 

128 m? 


(Poi ~Pr) 


(3) 


At steady-state, all the tube mass flows must add up to zero, and if all the tubes have 
the same diameter d and length /, the pressure differences must add up to zero as well. 

^rhi - 0 (4) 

1 = 1 


1 


( 5 ) 


X)(Poi-Pr) = 0 

1 = 1 

Substituting for p 0l in terms of the local velocity u u and multiplying by the presumed- 
constant tube spacing A y, we get 



EM - 

i = 1 

ipu, + iv. - !•’ ) &» = 0 

(6) 



“ ' m ') = 2(p ”“ - Pr) L 

1=1 

(7) 

where 


N 




L = ^2 Ay = N Ay 

(8) 


»=i 


is simply the height of the JV-tube rake as shown in Figure 1. 

The lefthand side summation in equation (7) is seen to be a midpoint-rule integration for 
the momentum+mass defect 

pul(t + S') 5 J(pu] - pu 2 w)dy (9) 

- - rf) &1K < 10 > 

l 

so that with the constant tube spacing Ay, the final result is 

pu 2 (0 + 8*) ~ 2 ( Pooo -Pr)L (H) 

or 0 + 8* ~ Po °° L (12) 

p 000 - Pe 

For the most accurate measurement of 9 + 8% it is clearly best to reference the reservoir and 
local-static pressure transducers to the freestream total pressure p 0oo , so that p 0oo ~Pr and 
p o —p e are direct measurements with no bias uncertainties. 

2 Near- wake Corrections 

Note that the integrating rake does not give the momentum thickness m isolation. If the 
measurements are to be used to support or validate drag calculations, this causes few prob- 
lems, since the measured sum 0 + 8* can simply be compared with the corresponding sum 
from the calculation. If the goal is absolute drag measurement, then it will be necessary to 
estimate the shape parameter H = 8* /0 to deduce the isolated 0. This can be done fairly 
accurately if only the minimum velocity in the wake is known. 

t/min = min ( 13 ) 
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Measuring this will typically require a separate total pressure tube placed sufficiently close 
to the wake centerline. 

Assuming the minimum velocity is known, a wake velocity profile is then quite closely 
approximated by the Coles cosine profile 


ny) 

u e 


Um'm "I" (1 Um'm ) 


1 ' 
2 . 


1 — cos 



—8 < y < 8 


where 8 is the wake half- thickness. From the definitions of 8* and 9, we have 


(14) 


8* = t ( l--)dy = (1 -11^)8 

J-t \ u e / 

9 = [‘ (--^)d y = (i-u^s -hi-u^fs 

J-S \u e U 2 e ) 4 

1 

a = 1 - 


(15) 

(16) 
(17) 


and the isolated momentum thickness then follows immediately from the measured B + 


8 * . 


9 + 8 
H + l 


The above results also apply to moderately assymetric wakes (such as the one sketched in 
Figure 1) since each wake half above and below the velocity minimum is still accurately 
represented by the Coles profile. 

Both U min and H quickly approach unity downstream, making their measurement or 
estimation less and less critical. Ideally, the rake should be moved downstream until it just 
captures the entire wake, but in practice the wake thickness and location are somewhat 
uncertain, and some conservatism will be required. 

However the shape parameter is estimated, a correction will still need to be applied to 
the resulting momentum thickness if u e /u „ at the rake is not unity, as is usually the case. 
The true profile drag/span is simply the ultimate momentum defect 

D' = pulB^ (19) 


where 9 ^ is the momentum thickness very far downstream. The evolution of the momentum 
thickness in the wake is given by the von Karman integral momentum equation 


i de 
9 dx 


~(H + 2 ) 


1 du € 
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( 20 ) 


which can be integrated from the rake position to far downstream if a unique function H(u e ) 
is assumed. The particular empirical assumption 


lnfito/ne) 
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V B-l J rake 


constant 


(21) 


results in the well-known Squire- Young formula. 
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( 22 ) 
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3 Design Requirements 

In light of the assumptions made in the derivation of relation (12), the following requirements 
must be met for the result to be valid. 

1. The flow in the tubes is fully- developed Poiseuille flow. 

2. The flow in each tube has negligible dynamic pressure. 

3. The static pressure at the tip of each tube is the local total pressure p 0 (y) ■ 

4. The pressure in the reservoir is uniform. 

5. The tubes all have the same length t, diameter d, and spacing Ay. 

Requirement 1 is equivalent to stating that the profile-development entrance length for 
each tube must be very small compared to the tube’s length. This gives the requirement 


or 


i > 


rhi 

47T/X 


l > 


d pu^d 

16 p 


(23) 

(24) 


which makes the worst-case assumption that p 0oo — p r — corresponding to the entire 

rake being immersed in a massive separated wake. In any case, it should be possible to meet 
the Poiseuille-flow requirement without too much difficulty by using sufficiently fine tubes. 
The only drawbacks of excessively-fine tubes are fragility, and possibly an excessive settling 
time. 

Requirements 2 and 3 are essentially equivalent, since negligible dynamic pressure in the 
tube automatically implies that the static pressure at the tip is nearly the full total pressure. 
The requirement is 



(25) 

(26) 


which turns out to be essentially the same as Requirement 1. 

Requirement 4 means that the dynamic pressures in the reservoir must be negligible 
compared to the pressure drops across the tubes. With the dynamic pressures in the tubes 
already being negligible, this additional requirement is met by making the cross-sectional 
area of the reservoir larger than all the tube areas put together. With D being the reservoir 
diameter, the requirement is 

D 2 > Nd 2 (27) 

which can be very easily met in practice. On the other hand, an excessively-large reservoir 
is undesirable, since it will increase the settling time. 

A design parameter which still must be selected is the appropriate overall height L. 
Clearly, this needs to be sufficiently large to capture the entire wake. However, making 
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it much larger than the wake width is undesirable, since increasing L will proportionately 
reduce the magnitude of the pressure signal p 0c0 — Pr, giving reduced accuracy. It is perfectly 
acceptable to build an oversize rake, and cap any tubes which fall well above or below the 
wake in order to improve accuracy. 

The remaining design parameter is the total number of tubes N . This only influences 
the accuracy of the midpoint summation approximation to the transverse profile integral. 
Twenty tubes should be sufficient for typical featureless wake profiles. If the wake thickness is 
expected to be significantly narrower than the rake height L, more tubes may be appropriate. 


4 Compressibility Effects 

If quantitative results are to be obtained from the integrating rake, corrections to the above 
relations may be necessary if the freestream Mach number is appreciable. The exact relation 
for the pitot total pressure is given by any of the standard isentropic formulas, e.g. 


Po(y) = Pe [^ + 1 Y L M 2 (y)Y 1 - Pei 1 


•» 2 (y) N 

2h a (y) / 


7-1 


( 28 ) 


For adiabatic flows at near-unity Prandtl numbers, the stagnation enthalpy can be assumed 
to be nearly constant across the wake and equal to its freestream value. 


h a (y) ~ h 0oo = 
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■1 Ml 
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Combining this with the second isentropic formula above gives 


Po(V) = Po a 


l + ^Ml 


u 2 (y)' 
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( 30 ) 


which is the relation which would be used with a conventional rake to determine the velocity 
profile from the measured total-pressure profile and the local static pressure p e . 

Just outside the wake, where u(y) = u e , the total pressure is the same as the freestream 
total pressure, which implies 


Pe_ 

Poo 


1 + 1 


u l). 


( 31 ) 


This gives an alternative form for the total pressure profile 


po(y) = p 0o 
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( 32 ) 
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which then replaces the equivalent incompressible form (2). A second-order Taylor series 
for equation (32) about = 0 conveniently recovers the incompressible form, but with a 
first-order Mach number correction. 
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1 - M 
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Po(V) p 0oo ~ lp~(u 2 e ~u 2 (y)) 
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2 
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+ 0(M 6 J (33) 


(34) 


Equation (34) truncates the Taylor series, and also drops terms which are O [( u\ - u 2 ) 2 } , 
which is justifiable for a typical small-defect wake profile. 

An additional effect of compressibility will be to produce a non-uniform kinematic viscos- 
ity profile across the wake. This will affect the tube mass flows Thj, which are still given by 
relation (3), but u must now be replaced by the local kinematic stagnation viscosity v 0 (y). 
Assuming that the dynamic viscosity varies as p ~ T b ~ h b results in 
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which is suitably approximated by its Taylor series. 
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T s = 110” K 


The viscosity exponent b can be estimated from the accurate Sutherland’s viscosity law 

, T . _ „ ( T T. + T s 
t>(T) - T + Ts 

by matching its logarithmic derivative at the freestream condition. 

3 1 


, _ dp 

H dT 


2 1 + Ts/T „ 


The rake tube mass flow relation (3) now becomes 

7T d 4 v. 


mi 


128 v e t u 0 i 

and the zero net mass flow condition (5) is 

(Poi-Pr) = 0 

i=i V « 


(. Poi - pr) 


or 
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(38) 
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The summation on the righthand side produces 


i # [fe - 0 L ~ ^ + «>] 


i=l 


where and are the kinematic thicknesses defined with the density profile omitted. 


9k + s * ~ /( 1_ ^) dy 


Substituting the total pressure approximation (34) into equation (42) produces 


^2—pJul-v^)^ 1 + ^Ml -f-1 Ay = 2(p 0oo -p r )E^ A y ( 45 ) 

^oi Poo _ V^oo /. Z = 1 ° l 


which then reduces to the compressible equivalent of equation (12) 
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where the higher powers of and ( u 2 e — u 2 ) 2 have been neglected as before. An explicit if 
somewhat imposing expression for the kinematic momentum+displacement thickness then 
follows. 
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Again, the natural reference pressure for the reservoir transducer is clearly p Qoo . The ratio 
of static/dynamic pressure in the equation above can of course be replaced by the Mach 
number 

-2=- = — (48) 

P.< -,Ml 

which is itself obtained from the total/static pressure ratio. The local edge/freestream 
velocity ratio follows from equation (31) 
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asked our local turbulence- measurement guru 
snny Breuer about the required rrdc sensitivity, 

.nee I wouldn't even hazard a guess. He's in 
reden , but I expect an email answer any minute now 



Appendix B: Key figures and sketches 



APEX 16 

x y 

1.000000 0.002036 

0.996936 0.002939 

0.988795 0.005338 

0.977291 0.008683 

0.962745 0.012817 

0.946684 0.017254 

0.930138 0.021685 

0.913447 0.026011 

0.896699 0.030208 

0.879909 0.034271 

0.863091 0.038196 

0.846244 0.041983 

0.829370 0.045629 

0.812468 0.049137 

0.795546 0.052501 

0.778591 0.055724 

0.761610 0.058805 

0.744595 0.061747 

0.727557 0.064550 

0.710498 0.067212 

0.693423 0.069734 

0.676338 0.072110 

0.659238 0.074342 

0.642136 0.076424 

0.625022 0.078355 

0.607898 0.080136 

0.590773 0.081760 

0.573635 0.083229 

0.556493 0.084542 

0.539344 0.085695 

0.522184 0.086690 

0.505023 0.087527 

0.487855 0.088202 

0.470681 0.088718 

0.453510 0.089073 

0.436337 0.089267 

0.419174 0.089299 

0.402019 0.089166 

0.384874 0.088867 

0.367749 0.088400 

0.350641 0.087762 

0.333560 0.086951 

0.316510 0.085961 

0.299492 0.084791 

0.282523 0.083434 

0.265599 0.081885 

0.248739 0.080140 

0.231948 0.078187 

0.215238 0.076023 


x y 

0.198625 0.073635 

0.182124 0.071014 

0.165761 0.068148 

0.149557 0.065024 

0.133556 0.061628 

0.117799 0.057944 

0.102355 0.053956 

0.087310 0.049650 

0.072785 0.045017 

0.058990 0.040088 

0.046262 0.034955 

0.035096 0.029839 

0.025944 0.025057 

0.018879 0.020836 

0.013576 0.017218 

0.009594 0.014109 

0.006567 0.011391 

0.004250 0.008957 

0.002489 0.006718 

0.001204 0.004607 

0.000381 0.002575 

0.000024 0.000624 

0.000089 -0.001269 
0.000556 -0.003190 
0.001466 -0.005090 
0.002906 -0.006846 
0.004935 -0.008327 
0.007560 -0.009522 
0.010851 -0.010544 
0.015033 -0.011514 
0.020527 -0.012515 
0.028026 -0.013603 
0.038406 -0.014821 
0.052060 -0.016156 
0.068128 -0.017542 
0.085370 -0.018957 
0.102840 -0.020404 
0.120161 -0.021923 
0.137275 -0.023559 
0.154200 -0.025357 
0.170985 -0.027355 
0.187738 -0.029589 
0.204639 -0.032060 
0.221631 -0.034644 
0.238067 -0.036998 
0.254021 -0.038916 
0.269935 -0.040388 
0.286070 -0.041481 
0.302451 -0.042255 


x y 

0.318997 -0.042742 
0.335659 -0.042958 
0.352439 -0.042918 
0.369357 -0.042645 
0.386421 -0.042174 
0.403596 -0.041531 
0.420860 -0.040737 
0.438200 -0.039812 
0.455601 -0.038770 
0.473059 -0.037626 
0.490567 -0.036393 
0.508116 -0.035082 
0.525706 -0.033703 
0.543327 -0.032267 
0.560976 -0.030780 
0.578654 -0.029249 
0.596354 -0.027682 
0.614079 -0.026084 
0.631831 -0.024461 
0.649604 -0.022820 
0.667405 -0.021167 
0.685216 -0.019511 
0.703005 -0.017861 
0.720758 -0.016229 
0.738477 -0.014622 
0.756158 -0.013049 
0.773800 -0.011521 
0.791402 -0.010046 
0.808960 -0.008636 
0.826477 -0.007301 
0.843943 -0.006051 
0.861354 -0.004901 
0.878696 -0.003863 
0.895951 -0.002958 
0.913086 -0.002205 
0.930043 -0.001634 
0.946688 -0.001280 
0.962660 -0.001176 
0.977023 -0.001319 
0.988508 -0.001619 
0.996818 -0.001946 
0.999998 -0.002103 





APEX 16 RE = 200000 MR = 0.600 Ncrit = 12.000 

APEX 16 RE = 200000 MR = 0.650 Ncrit = 12.000 

APEX 16 RE = 200000 MR = 0.700 Ncrit = 12.000 









APEX 16 RE = 500000 MR = 0.600 Ncrit =12.000 

RPEX 16 RE = 500000 MR = 0.650 Ncrit = 12.000 

APEX 16 RE = 500000 MR = 0.700 Ncrit = 12.000 






2000000 MR = 0.400 Ncrit = 12.000 

1000000 MR = 0.400 Ncrit = 12.000 











500000 MR = 0.650 Ncrit = 12.000 

300000 MR = 0.650 Ncrit = 12.000 








APEX 16 RE = 500000 MR = 0.700 Ncrit = 12.000 

APEX 16 RE = 300000 MR = 0.700 Ncrit = 12.000 













CDNT 



BL VELOCITY PROFILES 



